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Abstract. Radiation emitted from an accretion disk 
around a massive black hole is a widely discussed model 
for the UV/soft X-ray excess emission observed in the 
spectra of many AGN. A self-consistent calculation of the 
structure and the emerging spectrum of geometrically thin 
a-accretion disks in AGN is presented. The central object 
is assumed to be a Kerr black hole; full relativistic correc- 
tions are included. The local dissipation of gravitational 
energy is assumed to be entirely due to turbulence. Since 
these disks are mainly supported by radiation pressure, 
the assumption that the viscous energy production is pro- 
portional to the total pressure leads to diverging temper- 
ature structures in the upper parts of the disk, where the 
total optical depths is small. We therefore modify the stan- 
dard expression for the turbulent viscosity by taking into 
account the radiative energy loss of turbulent elements in 
an optically thin regime. Compton scattering is treated in 
the Fokker-Planck approximation using the Kompaneets 
oparator. The absorption cross section contains only free- 
free processes for a pure hydrogen atmosphere. We present 
several calculations for various parameters such as the ac- 
cretion rate M, the viscosity parameter a, the specific 
angular momentum a of the black hole and the inclina- 
tion angle ©o of the observer. The resulting temperature 
and density profiles show that the disks are optically thick 
for Compton scattering and effectively optically thin for 
most frequencies. This leads to spectra that are diluted 
with respect to the photon number but contain a Comp- 
tonizcd high energy tail in order to carry the required 
energy flux. In addition the electron temperature deviates 
strongly from the equilibrium temperature. For a model 
with M = 10 8 M Q , M = 0.3, a = 1/3 and a/M = 0.998 
the fraction of total flux emitted in the soft and hard X- 
ray band (> 2.4 • 10 16 iJz) for 9 = 0°, 41°, 60°, 70° and 
90° is 36%, 50%, 64%, 78% and 93%, respectively. There- 
fore, the model can in general account for the observed 
soft X-ray excess. 
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1. Introduction 

Dissipation of gravitational energy by accretion onto a 
massive black hole is regarded as the origin of the enor- 
mous luminosities of active galactic nuclei (AGN). Since 
the accreted matter is likely to carry some angular mo- 
mentum, it is common to invoke the presence of an accre- 
tion disk. A strong observational evidence for such disks 
comes from the spectral flattening in the optical/UV seen 
in many AGN, the so-called "big blue bump" (e. g. Shields 



1978; Malkan & Sargent 1982). This spectral feature is 



thought to be due to thermal emission from an accretion 
disk and perhaps extends into the soft X-ray band, pro- 
ducing a steep excess emission known as the "soft excess" 
(e. g. Arn aud e t al. |1985| ; Wilkes & Elvis |1987| Turner 



& Pounds 198E). Other evidence comes from the reflec- 
tion hump and iron emission line observed in AGN (e. g. 



Pounds et al. 1990) 



The standard theory of thin accretion disks is largely 
based on the fundamental papers of Shakura & Sunyaev 
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([1973J, hereafter SS73) and Novikov & Thorne ( |1973[ here- 
after NT73). The simplest way to calculate the emitted 
disk spectrum is to assume that the disk is geometrically 
thin but optically thick and radiates locally as a black- 
body. The effective temperature T e g of the blackbody is 
then solely determined by the dissipated flux per unit area. 
These simple disk models, however, are not able to pro- 
duce soft X-rays as long as the y do not become super- 
Eddington (e. g. Bechtold et al. 1987 ). Real disk spectra 
will differ from the blackbody approximation: for suffi- 
ciently high accretion rates and viscosity parameters, the 
accretion disk eventually becomes effectively optically thin 
at small radii and the gas temperature strongly deviates 
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from the equilibrium temperature in such cases. Even in 
the optically thick case the local spectra differ from the 
blackbody, because the scattering opacity dominates over 
absorption in the soft X-ray regime and there exists a 
temperature gradient in the vertical direction. 

In the last few years several calculations of increas- 
ing accuracy have been performed to determine the struc- 
ture and emission spectrum of accretion disks around mas- 
sive black holes. Contrary to the sta ndard Newtonian disk 
model, NT73 and Page & Thome ( |l974D calculated the 
effects of general relativity on the disk structure. The 
propagation of photons from the disk sur face to a dis- 
tant observer was treated by Cunni ngham ( |l975 ). On the 
other hand SS7 3, Czerny & Elvis ( |1987|) , and Wandel & 
Pctrosian (1988) concentrated on a proper calculation of 
the radiative transfer by discussing the effects of Comp- 
tonization in a simple analytic manner. Laor & Netzer 
( |1989| ) and Laor et al. (199C) included relativistic effects 
as well as free-free and bound-free opacities in their nu- 
merical computations. Most calculations of model spec- 
tra to date, however, made use of an averaging over the 



vertical direction. Ross, Fabian & Mineshige (1992, here- 
after REM) have calculated the vertical temperature pro- 
file and atomic level populations in the radiation pressure 
dominated inner regions of the disk for a given constant 
vertical density profile, using the Kompaneets equation 
to treat Compton scattering. Relativistic corrections for 
a non-rotating bl ack h ole were incorporated in this code 
by Yamada et al. ( 1994). Shi mura & Takahara ( 1993 ) and 
Shimura & Takahara ( 1995 ) have calculated the vertical 
structure and radiation field of a Newtonian disk self- 
consistently, using the ad hoc assumption that the local 
viscous heating rate is proportional to the mass density p. 

Here we have also calculated the z-structure simulta- 
neously with the radiation field of the disk. But we use a 
different viscosity description. We assume that the local 
energy production is caused by turbulence. The standard 
a-description (viscosity proportional to the total pres- 
sure) leads to diverging temperature profiles in the upper 
optically thin regions of the disk, because viscous heat- 
ing always exceeds radiative cooling. We therefore propose 
a modification of the standard a-model, which includes 
the radiative cooling of the turbulent elements. The fre- 
quency dependent radiative transfer equation is solved in 
the Eddington approximation and the effects of Compton 
scattering are treated by the Kompaneets equation (Kom- 
paneets |1957 ). Relativistic corrections on the local disk 
structure are introduced according to Riffert & Herold 
( |1995| ) for rotating and non-rotating black holes. In this 
paper we present solutions of the local structure and the 
emission spectrum of accretion disks for different input pa- 
rameters (accretion rate, viscosity parameter). Integrated 
disk spectra, as seen from a distant observer, are calcu- 
lated by the use of a transfer function (e. g. Cunningham 
1975 ; Speith et al. 1995 ). In Sect. 2 the basic equations 
and viscosity description of our model are formulated. In 



Sect. 3 numerical results are presented. Finally, our results 
and conclusions are summarized in Sect. 4. 



2. The Model 

2.1. Basic Equations 

We adopt the standard geometrically thin a-accretion 
disk model, i.e. the disk height H is much smaller than 
the radius r, and the dominant velocity is given by the 
Keplerian motion around the central mass. The disk is as- 
sumed to be in a stationary and rotationally symmetric 
state, thus all functions depend only on the radial and 
vertical coordinates r and z. The relativistic disk struc- 
ture has been calculated by NT73 and subsequently by 



Riffert & Herold (1995), correcting a term in the vertical 
pressure balance. According to this paper we define four 
relativistic correction factors with respect to the standard 
Newtonian disk model depending on the mass M and the 
specific angular momentum a/M of the central black hole: 
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Ti is the inner disk radius located at the position of the 
last stable circular orbit; the gravitational constant and 
the velocity of light are denoted by G and c. 

The hydrostatic equilibrium in vertical direction for a 
thin accretion disk is given by 



dP 

dz 



-P9z 



-?- 



GM C 



.3 Z B > 



(2) 



where P, p, and g z are the total pressure, the mass density, 
and the gravitational acceleration, respectively. 

The radiative transfer is solved in the Eddington ap- 
proximation, and Compton scattering is treated in the 
Fokker-Planck approximation using the Kompaneets op- 
erator. The absorption cross section contains only free- 
free processes for a pure hydrogen atmosphere. Induced 
contributions to the radiative processes have been ne- 
glected throughout. 

Assuming the plasma to be in a state of local thermo- 
dynamic equilibrium (LTE) the first two moments of the 
radiative transfer equation for the spectral energy density 
U v and the spectral flux F v read 



dF v B 8irh 2 4 
- K u pc[W v -U v ]+K T p TV— v 



dz 



A n 



kT dn u 



h dv 



(3) 
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1 BU V 
3 dz 



(4) 



Here W v = 87r/ii/ 3 c -3 exp(— hv/kT) is the Wien function 
which serves as the equilibrium spectrum in this case, ^ 
is the free-free opacity, and kt = 0.4 g/cm 2 is the Thom- 
son opacity. The photon occupation number n v can be 
expressed in terms of U v 



The free-free opacity is given by 
= 1.3 • 10 56 pT- 1/2 v- z g s cm 2 g' 1 



(5) 



(6) 



where T is the gas temperature, k is the Boltzmann con- 
stant, and <?ff is the thermal average free-free Gaunt factor 



(Karzas & Latter 1961; Carson 1 



In the Eddington approximation the equation of state 
can be written as 



P = Pn, 



Pn 



pm u 



-pT 



1 



(7) 



where U = J °° U v dv is the total energy density of the 
radiation field, m u is the atomic mass unit and \i is the 
mean molecular weight with fj, = 1/2 for a fully ionized 
hydrogen atmosphere. 

For a Keplerian disk, the viscous heating rate per unit 
volume Q V i S c is related to the </>-r-component of the vis- 
cous stress tensor t by 



3 IGM A 



(8) 



and t^r can be expressed in terms of velocity gradients 
leading to the result 



3 GM A 

2 v y—p 



(9) 



with the shear viscosity rj. 

In a steady state this heating rate must be balanced 
by the total radiative cooling, since convection is not con- 
sidered in this model. Integration of the right-hand side of 
equation (||) over frequency v leads to the following energy 
balance equation 



dF 

dz 



dz 



dv . 



(10) 



The disk model is so far described by the equations (0), 
@) @> © and © which can be solved numerically 
when proper boundary conditions are imposed. The only 
unspecified process is the source of viscosity and this will 
be discussed in the next section. 

Because of symmetry reasons the inner boundary con- 
dition in the midplanc of the disk (z = 0) is taken as 



F v {z = 0) = , 



(11) 



and at the upper boundary (at z = H) we assume that 
the disk radiates isotropically into the vacuum 



F v = -V v . 



(12) 



For a geometrically thin accretion disk the energy flux 
from the disk surface (z — H) is 



/•OO 

F(z = H)= F v dv 
Jo 



3GMM D 
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(13) 



where M is the total mass flux through the disk. 

Finally, the disk height H is specified by the condition 



p{z = H) = 



(14) 



We solve the above set of equations for the functions p, 
T, H, U„, and F v using a finite difference scheme in both 
variables z and v. The vertical structure is resolved with 
100 points on a logarithmic grid, and 64 grid points are 
used in frequency space. The resulting set of algebraic dif- 
ference equations is then solved by a Newton-Raphson 
method. 

Once the vertical structure has been determined we 
can integrate the density profiles to obtain the surface 
density £(r) = p{z) dz and the average radial veloc- 
ity \v r \ — M/S(r). In order to justify the basic assump- 
tions of the thin disk model one must have \v r \ <g; \v$\, 
where v^, is the Keplerian velocity. 

2.2. The Viscosity Description 

In the standard a-model (SS73) the viscous stress tensor 
component tsr is propotional to the pressure 



tfr. = aP 



(15) 



This is a consequence of a number of assumptions: first, 
the viscosity is caused by turbulence in the disk, and the 
magnitude of 77 can be estimated from a dimensional anal- 
ysis 



V ~ P hurbVturb 



(16) 



where hurb and Vturb are the size and the typical veloc- 
ity of the largest turbulent eddies, respectively. Next, the 
disk height H is used for Uurb which is obtained from the 
hydrostatic equation (||) 



(17) 



where P c is the pressure in the midplane, and p stands 
for the average density. The turbulent velocity Vturb is 
assumed to be limited by the local sound speed c s . Oth- 
erwise shocks would develop and heat the gas until the 
tubulcncc is subsonic. Approximating the ratio P c j p by 
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the local value Pj p — c 2 s and inserting the above relations 
into equation (^|) for the viscous shear leads to the expres- 
sion ([l5]), where the constant a takes up all uncertain- 
ties in terms of numerical factors. This entire procedure 
is quite successful for the standard disks in cataclysmic 
variables which are dominated by the gas pressure. If ra- 
diation pressure becomes important a number of difficul- 
ties arise when the model is applied without modifications. 
According to equation (||) the radiative cooling decreases 
in the upper layers of the disk where the density drops ex- 
ponentially. Since the radiation pressure remains finite in 
this region, the total heating rate (@) cannot be balanced 
by radiative losses. This could be cured by replacing P in 
equation ( |l5| ) by the gas pressure, however, in an optically 
thick regime the radiation behaves like a gas exchanging 
momentum with the plasma and thus contributes to the 
viscous shear. In addition, in a radiation dominated disk 
gravity is mainly balanced by radiative forces and the disk 
height is then given by 



H 



3Mk f D 



ire C ' 

where kf is a flux weighted average opacity 



1 

Hp 



(18) 



(19) 



This height does not contain the sound speed. In order 
to derive a consistent parameterization of the turbulent 
viscosity for the radiation dominated case we start from 
the dimensional analysis (^|) using again Iturb = H, where 
the disk height H is self-consistently calculated in our 
model. For the limitation of the velocity Vturb the same 
arguments as in SS73 are used, i.e. Vturb has to be small 
enough that no entropy due to shock waves is generated in 
the flow. For a radiation dominated shock in an optically 
thick medium 



with the optical depth 



/ Kppdz' . 

J z 



(23) 



and the pressure ratio (3 = P gas /P- For the turbulent 
viscosity we now have 



77 = apHVt 



t uro 



(24) 



with a constant a which is a parameter in our disk model. 

We will now consider the heating-cooling balance, in- 
tegrate equation (||) over frequency, and use the relations 
m and (|) 



dF 



R s pc [a r T 4 -U]+4U — [kT- <v>] 
dz m e c 



9GM A 2 



apHVt 



turb 



(25) 



where Kg is some appropriately averaged mean absorption 
opacity, a r is the radiation constant, and 



<v>= 



/ °° u A n v dv 
4 J °° v 3 n v dv 



(26) 



is an average photon energy (h < v > / 'k = T r can be 
considered as a radiation temperature). In the upper parts 
of the disk (at z w H) the density drops to zero and since 
Kff oc p free-free processes are negligible in this regime. 
The radiation field remains unchanged and thus < v > and 
U are constant; from the boundary condition cU(H) = 
2F(H) we obtain using equations (|l3|), ([l8|), and f^ ) for 
t < 1 



T — T r = , 



(27) 



Vturb < 



(20) 



but in an optically thin regime, where the photons are able 
to escape freely 



Vturb < 



Pqas 



P 



(21) 



This follows from considering an isothermal shock in a 
plasma with radiation when the Mach number approaches 
unity. Assuming an isothermal shock structure is an ap- 
proximation for an optically thin environment. Since this 
limiting value of Vturb depends on the optical depth we 
interpolate between the two extreme cases 



*turb C s 



VP 



l + T 



where 



r -'i)'(^' 



2 2 

ra„<r 



pkm v 



= 5.1 10 b a 2 



^) K.(28) 



(Note that in a scattering dominated plasma kf/kt ~ 1). 
From this we get the temperature that will result from a 
balance of viscous heating and inverse Compton cooling 



T = T r + \ 



T v + y/T* + 4T r T v 



(29) 



This solution is unique since a negative sign of the square 
root in ( p9|) would correspond to a negative sign of the 
right hand side of equation (^7|). As a consequence we 
obtain a lower limit for T 



(22) T > max[T r , T v ] 



(30) 
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3. Results 

In this paper the mass of the black hole is fixed at a typ- 
ical value of M = 10 8 M Q . We show results of our model 
calculations for different input parameters such as the ac- 
cretion rate M, the viscosity parameter a, the specific 
angular momentum a, and the inclination angle Qq of the 
observer's position with respect to the disk axis. Here- 
after, M is measured in units of the critical accretion rate 
M cr it = L Edd /(ec 2 ), where L E dd = 4ttcGM/k t is the Ed- 
dington luminosity and e is the efficiency of accretion. For 
a non-rotating black hole (a/M — 0) e is 0.057, and for a 
maximally rotating black hole (a/M ~ 0.998, see Thorne 
1974) e ~ 0.321. First, we give solutions of the local ver- 
tical structure and emission spectrum at 5 Schwarzschild 
radii (r — 5Rs — lQGM/c 2 ). At this radius, the local en- 
ergy release takes its maximum value for a non-rotating 
black hole and this disk region is responsible for the high 
frequency part of the total spectrum. Next, we consider 
the radial structure and the integrated disk spectrum seen 
by a distant observer. At this point we also include rela- 
tivistic effects on the emergent disk spectrum such as the 
relativistic Doppler shift from the disk rotation, gravita- 
tional redshift, and gravitational light bending due to the 
central black hole. 



3.1. Vertical Structure and Local Disk Spectrum 

A typical vertical distribution (normalized to the corre- 
sponding maximum values) of the radiative flux F, the 
radiation energy density u, the mass density p, and the gas 
temperature T as a function of the Thomson scattering op- 
tical depth tt is shown in Fig. |l| for a non-rotating black 
hole with M = 1O 8 M , M = 0.3, a = 1/3 and r = 5R S . 
In Fig. |2|, the above variables are plotted for the case of 
a maximally rotating black hole (a/M = 0.998) with the 
same input parameters, but using r = IRs , which nearly 
corresponds to the maximum of the local energy release 
for a/M = 0.998. It is found that the mass density p is 
not a monotonically decreasing function of z, but is show- 
ing a clear density inversion (such density invers ions h ave 
also been fou nd by Meyer & Meyer-Hofmeister 1982 and 
Milsom et al. 1994 ). This can be explained as follows: in 
using the flux weighted mean Kp (equation (|l9|)), the hy- 
drostatic equilibrium can be written as 
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Fig. 1. Vertical distribution of radiative flux F, radiation en- 
ergy density U, mass density p, and gas temperature T (all 
normalized to their respective maximum values) at r = 5Rs- 
The model parameters are: M = 10 8 M Q , M = 0.3, a = 1/3, 
a/M = 0. Note that U approaches a finite value for small op- 
tical depth (tt — > 0). 




100.00 10.00 



Fig. 2. Same as Fig. 1, but for a maximally rotating black hole 
a/M = 0.998 at r = 1R S . 



dP n , 



dz 



p[(F(z) 



H Edd 



(*)] 



~P9z 



where the local Eddington flux FEdd(z) is given by 
c 



F E dd(z) 



up 



(31) 



(32) 



and g* = g z — KpF/c is the effective gravitational ac- 
celeration. For a highly radiation pressure dominated disk 
(Prad ^S> Pgas), the local flux adjusts itself close to the Ed- 
dington value. For example, the radiation pressure always 



exceeds the gas pressure by at least a factor of about 100 
throughout the disk for the parameters adopted in Fig. 
0. When the emerging surface flux is large, the local flux 
becomes slightly super-Eddington in the inner regions of 
the disk and a strong density inversion (dp/dz > 0) occurs. 
Because of the boundary condition imposed on the surface 
flux, the outer parts of the disk are always sub-Eddington. 
In Fig. |!| the vertical density distribution is shown for dif- 
ferent M. The density in the equatorial plane is up to a 
factor of 7 smaller than its maximum value. Although the 
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Fig. 3. Vertical density distribution at r = 5Rs for a model 
with a = 1/3 and a/M — 0. The curves are labeled with the 
value of M. 



total Thomson scattering depth decreases for increasing 
M and a, our model calculations show, that the accretion 
disks are optically thick with respect to Thomson scatter- 
ing even for large M and a. At the same time these disks 
are effectively optically thin {r e s{v) < 1) except for the 
low frequencies, and T c g(v) is defined as 



H 



T e s(v) = / /V 3 ■ [«!ff (y) + k t ]k s (v) dz . 
Jo 



(33) 




In this case the assumption of local thermodynamic equi- 
librium (LTE) is no longer valid even in the disk midplane 
and the gas temperature deviates strongly from the equi- 
librium temperature defined as T eq — (U/dr) 1 / 4 . The var- 
ious heating and cooling functions are compared in Fig. 
f|. Compton cooling plays the dominat role in the inner 
layers of the disk (at small z), whereas in the region of the 
density increase free-free-processes become very efficient 
leading to a steep temperature drop. The gas temperature 
even falls below the radiation temperature T r , thus Comp- 
ton scattering turns into a heating mechanism (A comp is 
negative) which effectively balances the bremsstrahlung ^ 
cooling. In the uppermost layers, where the density drops 
exponentially, free- free processes are negligible (oc p 2 ), and 
the temperature adjusts itself to the equilibrium between 
viscous heating and inverse Compton cooling described 
above. Due to these heating-cooling processes a temper- 
ature inversion occurs and a hot corona-like layer builds 
up on top of the disk. This layer is however optically thin 
to both scattering and absorption (the typical densities 
are below 10 -13 gcm -3 ), and thus there is no effect on 
the emitted spectrum. Note that in the upper disk lay- 
ers the viscous heating is proportional to the density, thus 
the relevant cooling mechanisms must not have a stronger 
p-dependence. This holds of course for inverse Compton 



0.22 0.24 0.26 0.28 0.30 0.32 0.34 
z [RJ 

Fig. 4. Viscous heating A V i SC , Compton cooling A comp and 
cooling due to free-free emission Ag (all normalized to the 
viscous heating rate at the disk midplane) as a function of z at 
r = 5Rs- The model parameters are: M = 1O 8 M , M = 0.3, 
a — 1/3, a/M — 0. The outermost layer is shown in the upper 
right corner with higher resolution. 



cooling used in our model, however, cooling due to strong 
resonance lines has also been proposed (Hubeny 199C ) to 
contribute to the overall heating-cooling balance. 

The gas temperature and the corresponding equilib- 
rium temperature T eq for different M are shown in Fig. ^. 
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Fig. 5. Vertical distribution of gas temperature T (solid lines) 
and equilibrium temperature T eq (dotted lines) at r = 5Rs for 
a model with a = 1/3 and a/M = 0. The curves are labeled 
with the value of M. 
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For large values of M even the temperature in the mid- 
plane of the disk strongly exceeds the equilibrium value. 
For low M the disk becomes effectively optically thick and 
a general equilibrium is established. 

The locally emitted spectrum at a fixed radial position 
(r = 5i?s) for different M is shown in Fig. |^. The local 




Fig. 6. Emergent spectrum at r = 5Rs for a — 1/3 and 
a/M — 0. The curves are labeled with the value of M. Dotted 
curves show the local Wien spectrum at T = T e ff. 



Wien spectra at the corresponding effective temperature 
T e ff are also shown by dotted lines. At low frequencies, the 
accretion disk is effectively optically thick and scattering 
effects can be neglected; the local spectrum approaches 
the equilibrium distribution W u , however, the spectra in 
low frequency regime cannot be taken too seriously since 
we have neglected induced processes. At higher frequen- 
cies, scattering effects become important and one gets a 
modified Wien spectrum. For sufficiently high M and a 
a Wien peak appears in the soft X-ray range of the local 
spectra. This feature is a result of repeated Compton scat- 
tering of photons by a thermal distribution of electrons at 
temperature T. If the Compton parameter i/> 1, Comp- 
tonization goes to saturation and the photons are shifted 
into a Wie n distribution at temperature T (see Felten & 
Rees 1972). The Compton parameter y is defined as 



AkT 



y 



'T ■ 



(34) 



For example, the case M — 0.3 shown in Fig. ^ leads 
to a effective optical depth T e g(v max ) = 0.1 at the peak 
frequency v max = 2.8 • 10 16 Hz and a total Thomson scat- 
tering depth of tt = 98.5; from that and an average tem- 
perature of T = 4 • 1Q 5 K (sec Fig. ||) we have y » 3. Since 
repeated Compton scattering is important for y > 1, a sig- 
nificant shift of photons into the Wien peak is expected. 



For low M and a the disk is effectively optically thick 
for all relevant frequencies and the radiation assumes an 
equilibrium spectrum for almost all z. The effective op- 
tical depth for the curve M = 0.1 in Fig. || and |6| is 
T c s{v max ) = 9.8, with v max = 9.1 ■ 10 15 Hz. Neverthe- 
less, the emergent spectrum deviates from the Wien spec- 
trum at T — T ff because there exists a temperature gradi- 
ent and high frequency photons escape from deeper layers 
with higher temperatures. 

3.2. Radial Structure and Integrated Disk Spectrum 

After discussion of the vertical structure and local emis- 
sion spectrum at r = 5Rs here we investigate the radial 
structure and overall spectrum of the accretion disk. In 
numerical calculations, we used 50 radial points on a log- 
arithmic grid from the last stable orbit rj to the outer disk 
radius r out (here r out — 1000i?s). The height of the disk 
from the surface to the equatorial plane as a function of 
the radial distance from the black hole is shown in Fig. [?] 
for several values of M. For comparison, the geometrically 
thin limit, here defined as H(r) — O.lr, is also plotted as 
a dashed curve. This means that the accretion rate should 




R [Rs] 

Fig. 7. Disk height for a = 1/3 and a/M = 0. The curves 
are labeled with the value of M. The dashed curve shows the 
geometrically thin limit H(r) = 0.1 ■ r. 



not greatly exceed M = 0.3 for not to violate the geo- 
metrically thin disk approximation. Therefore, numerical 
results for higher accretion rates are strictly speaking not 
self-consistent. On the other hand, the a-parameter has 
little influence (of the order of a few per cent) on the 
height of the disk. If the local flux is given by the Ed- 
dington flux, the disk height H is given by equation ([l8|). 
In the inner regions of the disk Thomson scattering de- 
livers the major contribution to the opacity (kf ~ Kt) 
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and the height of the disk is solely determined by the ac- 
cretion rate M. However, the local flux deviates from the 
Eddington flux at the upper boundary of the disk. Al- 
though the gas pressure is negligible in comparison with 
the radiation pressure, the gas pressure gradient dpcas/dz 
is not. We have also determined mean radial velocities 
|u r | = M/S(r) for all radial grid points. Our calculations 
show, that for M < 0.3 and a < 1 the basic thin disk 
approximation |w r | <C \v^\ is always satisfied. Now we 
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10' J 10' 
V [Hz] 

Fig. 8. Entire Newtonian face on disk spectra for a — 1/3 and 
a/M = 0. The curves are labeled with the value of M. Dotted 
curves denote multi Wien spectra. 
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JO 




10' J 10' 
V [Hz] 

Fig. 9. Entire Newtonian face on disk spectra for M = 0.3 
and a/M — 0. The curves are labeled with the value of a. The 
dotted curve denotes the corresponding multi Wien spectrum. 



will turn to the integrated disk spectra. General relativis- 
tic effects on the emergent spectrum are twofold: first, 
relativistic effects occur in calculating the disk structure 
and therefore change the local emission spectrum in the 
corotating frame of the disk. The resulting spectra were 
discussed for a non-rotating black hole in section 3.1. Sec- 
ond, relativistic effects have a substantial influence on the 
propagation of photons to the observer. All these effects 
become more and more important if the black hole is ro- 
tating, since then the disk extends to regions very near to 
the black hole. To separate both effects the investigation 
of the integrated disk spectrum is treated in two different 
ways: first, we examine the emergent spectrum of a face on 
disk as seen for a distant observer without including rela- 
tivistic effects or Doppler shifts on the locally emitted pho- 
ton distribution. Relativistic effects only occur due to the 
relativistic disk structure equations (see section 2.1). The 
entire spectrum then is simply calculated by integration 
of the local spectra over the disk surface. In the following 
we will call this the Newtonian approximation. Second, a 
fully relativistic calculation of the propagation of photons 
from the disk to a distant observer is performed, using 



a program code (Speith et al. 1995) to obtain numerical 
values of the Cunningham transfer function (Cunningham 



1975) for any set of parameters. 



The Newtonian disk spectrum for several values of M 
and a is shown in Fig. |^ and ^. The spectral luminosi- 
ties calculated by summing up local Wien spectra at the 
corresponding effective temperatures are indicated by dot- 
ted lines. The hardness of the spectra and therefore also 
the amount of soft X-rays is a sensitive function of M 
and a somewhat less sensitive function also of a. In the 
following, we fix the accretion rate and viscosity parame- 
ter at M = 0.3 and a — 1/3, respectively. A comparison 
of the relativistic and Newtonian face-on disk spectrum 
for a non-rotating (a/M = 0) and maximally rotating 
(a/M = 0.998) black hole is shown in Fig. [l(]. For low 
frequencies the emergent spectrum approaches the New- 
tonian model. The spectrum at higher frequencies is re- 
duced with respect to the Newtonian case especially for 
a rotating black hole. These photons originate from the 
inner parts of the disk with high velocities and gravita- 
tional fields and therefore are diverted by forward peaking 
and gravitational focusing. Consequently, an equatorial 
observer primarily sees blueshifted and focused radiation 
from the hot inner parts of the disk. In Fig. [ll] the observed 
relativistic spectrum is shown as a function of the incli- 
nation angle Qq together with a Newtonian face-on spec- 
trum (dotted line) for a maximally rotating black hole. A 
significant fraction of total flux is emitted in the soft and 
hard X-ray range (v > 2.4- 10 16 Hz) and strongly depends 
on the position of the observer. For 9 = 0°, 41°, 60°, 70° 
and 90°, this fraction is 36%, 50%, 64%, 78% and 93%, 
respectively. The fraction of the total flux emitted in the 
soft X-ray band ( 2 . 4 • 1 16 iiz - 6 • 1 1 7 iiz ) , corresponding 
to the sensitivity range of the PSPC (position sensitive 
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Fig. 10. Entire face on disk spectra for M = 0.3 and a — 1/3. 
The curves are labeled with the value of the specific angular 
momentum a/M. The solid an dotted curves denote the rela- 
tivistic and Newtonian disk spectra, respectively. 



1.0 




Fig. 11. Entire relativistic disk spectra for different inclination 
angle 9 = 0°, 41°, 60°, 70° and 90°. The model parameters 
are: M = 10 s M , M = 0.3, a = 1/3, a/M = 0.998. The dotted 
curve denotes the Newtonian face-on spectrum. 



proportional counter) on the ROSAT X-ray satellite is 
32%, 38%, 38%, 33% and 21%, respectively, and therefore 
is nearly independent of the inclination angle. In order to 
show the limiting case of our calculations we also included 
cos 0o = in Fig. ffl[ However, for a nearly edge-on ob- 
server self-occultation of the inner parts of the disk by the 
outer parts should be taken into account. 



4. Concluding Remarks 

We have performed a self-consistent calculation of the ver- 
tical structure and emergent spectrum of an accretion disk 
around a massive Kerr black hole. Full relativistic correc- 
tions have been included. The standard a-model leads to 
diverging temperature profiles in the upper layers of the 
disk, where the density drops exponentially and viscous 
heating always overcomes radiative cooling. When we in- 
clude the radiative cooling of the turbulence elements in 
the optically thin part of the disk, an equilibrium between 
turbulent viscous heating and inverse Compton cooling is 
established. Our calculations are for a fixed central mass 
M = 1O 8 M but for different M, a, a/M, and 6 . It 
is found that the mass density is not a monotonically de- 
creasing function of z but shows a strong density inversion 
in the radiation pressure dominated inner parts of the disk, 
where the local flux becomes slightly super-Eddington. 
Moreover, we always obtain a temperature inversion in 
the upper optically thin layers of the disk with zero gradi- 
ent in the outermost region. For sufficiently high M and 
a the disk becomes effectively optically thin at small radii 
and the gas temperature exceeds the equilibrium temper- 
ature (T eq = (U/cir) 1 / 4 ) by a large amount even in the 
disk midplane. As a result, the local emission spectrum 
strongly deviates from the Wien spectrum at T = T e g. 
Even in the optically thick case, the local spectra dif- 
fer from the Wien-spectra at high frequencies, because 
scattering opacity dominates in the soft X-ray range and 
there exists a temperature gradient in the vertical direc- 
tion. Therefore, a significant fraction of the total flux is 
emitted in the soft X-ray range even for low M and a, 
especially in the maximally rotating case (a = 0.998) and 
for high inclination angles of the observer. For increasing 
M, a, a/M, and ©o the calculated spectra become more 
and more flat, producing quasi power law spectra in the 
sensitivity range of the ROSAT PSPC, whereas for small 
values a steep soft X-ray component is established. Thus, 
the model can in general account for the soft X-ray ex- 
cess observed in many AGN. We intend to show detailed 
spectral fits of this model to combined ROSAT and IUE 
observations for a sample of AGN in a future paper. 

Note that we have neglected induced processes and 
that our opacity description only contains free-free pro- 
cesses for a pure hydrogen atmosphere. Neglecting induced 
processes means, that we have concentrated on the high 
energy tail of the resulting spectra. Nevertheless, bound- 
free absorption and line opacities can significantly con- 
tribute to the total opacity especially at low temperatures. 
Therefore, the spectrum is expected to soften when all 
relevant opacities are taken into account because of the 
increasing effective optical depth in this case. 

Another important point is the role of convection in 
our calculations, especially in view of our inverse density 
profiles. If convective energy transport takes place in the 
disk, the vertical structure and subsequently the emergent 
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spectrum may be altered considerably. We therefore have 
done a simplified stability analysis according to the stan- 
dard mixing length theory (see for example Cox & Giuli 



1968), assuming that the convective elements move adia- 



batically in the disk and that the local gas temperature is 
given by the equilibrium temperature T eq . In a gas with 
radiation pressure, the effective buoyancy force is reduced 
by the radiation pressure gradient and the convective flux 
then can be expressed in terms of the effective gravita- 
tional acceleration g* (equation j3l|) , 



± conv — 



pc P TllH-^ 2 5^ 
4v^2 



v/<? z *(V s -V e )-(V s -V e ),(35) 



where l rn , Hp — P/(pg z ), and cp are the mixing length, 
the pressure scale height and the specific heat at con- 
stant pressure, respectively, and <5 = —dlnp/dlnT. V s = 
(dlnT / dlnP) s and V e = (dlnT / dlnP) e are the tempera- 
ture gradients with respect to pressure for the surround- 
ings and the rising element. Therefore, the disk is unsta- 
ble for convection if the square root on the right hand 
side of equation ( |35|) has a positive argument. In the sub- 
Eddington regime one has g* z > and the condition for 
convection to set in is given by the Schwarzschild criterion 
(V s — V e ) > 0. On the other hand in the super-Eddington 
regime g* < 0, and the condition for convective instabil- 
ity is (V s — V e ) < 0, leading to a convective energy flux 
in negative z-direction. Our model calculations show that 
none of both conditions are satisfied and the disk there- 
fore is stable against convection. Since the gas temper- 
ature deviates from the equilibrium temperature for an 
effectively optically thin disk and the convective elements 
do not move adiabatically in the upper optically thin lay- 
ers, this simplified stability analysis has to be considered 
with caution, and a more detailed stability analysis has to 
be performed in order to clarify the role of convection. 

It is well known that an a-disk with a viscosity propor- 
tional to the total pressure is radially unsta ble in the ra- 
diation dominated case (Lightman & Eardly 1974 ). How- 
ever, this analysis is based on a vertically integrated disk 
structure. Since our model leads to a strong z-dependence 
for almost all functions including the viscosity, a detailed 
2-dimensional stability analysis would be required for a 
definite answer concerning the overall stability of our 
model. 
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